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Abstract 

We investigate magnetism and quantum phase transitions in a one-dimensional system of in- 
tegrable spin-1 bosons with strongly repulsive density-density interaction and antiferromagnetic 
spin exchange interaction via the thermodynamic Bethe ansatz method. At zero temperature, 
the system exhibits three quantum phases: (i) a singlet phase of boson pairs when the external 
magnetic field H is less than the lower critical field Hd] (ii) a ferromagnetic phase of atoms in 
the hyperfine state \F = \,mp = V) when the external magnetic field exceeds the upper critical 
field Hc2] and (iii) a mixed phase of singlet pairs and unpaired atoms in the intermediate region 
Hci < H < Hc2- At finite temperatures, the spin fluctuations affect the thermodynamics of the 
model through coupling the spin bound states to the dressed energy for the unpaired mp = 1 
bosons. However, such spin dynamics is suppressed by a sufficiently strong external field at low 
temperatures. Thus the singlet pairs and unpaired bosons may form a two-component Luttinger 
liquid in the strong coupling regime. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.IK, 05.30.Fk 



I. INTRODUCTION 



The well-developed techniques for controlling and manipulating Bose-Einstein conden- 
sates (BECs) of spinor atoms provide an excellent opportunity to explore novel magnetism 
and quantum phases. In a magnetic trap, atoms with different magnetic moments are sub- 
jected to different forces, so it is very difficult to confine a spinor BEC involving all possible 
spin states. However, since the laser-atom interaction is determined by the induced electric 
dipole moment, an optical trap may confine every spin state to preserve the "vector" prop- 
erty of spinor atoms. This allows one to trap a true spinor BEC, which involves an ensemble 
of Bose atoms condensed in a coherent superposition of all possible hyperfine states. In this 
way, several experimental groups have successfully demonstrated spinor BECs of ^'^Na 
and ®^Rb js, 4| atoms. 

The ground states and some low-energy excitations of a spinor BEC were theoretically 
analyzed by Ohmi and Machida [sl and Ho j^. It has been shown experimentally that all 
three spin components of a spin-1 condensate can be either miscible or immiscible with one 
another where the immiscibility will lead to formation of spin domains This has been 
confirmed numerically using the Gross-Pitaevskii equation and Thomas-Fermi approxima- 
tion j^. Using the single-mode approximation (SMA), the ground state population dynamics 
of a spin-1 BEC have been studied by Law et al. [s] and Pu et al. [10]. They have also found 
that the ground state is a superposition of collective spin states (Fock states) and cannot 



be expressed clS 9j product of individua 
all three spin components. Ho and Yip 



_spin states. This shows the collective behavior of 



11| , found the antiferromagnetic ground state to be 



a fragmented condensate with large particle number fluctuations as stated in the references 







9|, |10|. This fragmented condensate gradually deforms into a more stable coherent state 
as the strength of the external fleld gradient increases. Recently, Rizzi et al. [12] applied 
the DMRG method to determine the phase diagram for one-dimensional spin-1 bosons. In 
accordance with Imambekov et al. 13] and Yip [l^, they showed that the dimerized state 
is among the ground states. The quantum phases such as the polar phase, nematic phase 
and spin singlet phase were discussed by Demler and Zhou [l5|. 



In one dimension (ID), spinor Bose gases have a ferromagnetic ground state in the ab- 



sence of spin-dependent forces 



181 ]. However, the spinor Bose gas can have either 



a ferromagnetic or an antiferromagnetic ground state in the presence of spin-exchange in- 
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teraction 



fly , ll9|]. Very recently, Cao et al. [20|] 



191 ] ■ Very recently, Cao et al. [20|] proved that there exists an integrable 



point in scattering parameter space for ID spin-1 bosons with both delta-function contact 
interaction and spin-exchange interaction. This model provides an important benchmark to 
understand spinor BECs and spin liquids in low dimensions. From the exact Bethe ansatz 
(BA) solution, Cao et al. found that the ground state is a spin singlet in the absence of 



an external field. Essler et al. 



21| then proved that the low energy physics in the weak 



repulsive coupling regime can be described by a spin-charge separated theory of an effec- 
tive Tomonaga-Luttinger Hamiltonian and an 0(3) nonlinear sigma model. In this weak 
coupling limit, both the collective pairing fluctuations and spin fluctuations dominate the 
low-lying excitations. Using the BA equations and the effective field theory, Essler et al. 
calculated the scaling dimensions and the large-distance asymptotics of correlation functions 
of the model. Such a spin liquid phase was also previously investigated by Zhou |22| through 
the introduction of the Weyl representation of SU{2). 

In this paper, we investigate quantum liquid phases in the ID integrable system of spin-1 
bosons 20|] with strongly repulsive and antiferromagnetic spin-exchange interactions. We 
derive the thermodynamic Bethe ansatz (TBA) equations on the basis of particle-hole exci- 
tations 3] and spin-strings 24 1. From these equations we find that for the strong coupling 
regime there is a large energy gap in the lowest spin excitation. We also show that spin 
fluctuations are frozen out under a strong external field at zero temperature. When the 
external magnetic field exceeds the lower critical field ifci, the energy gap vanishes and the 
charge excitations evolve into two gapless modes of singlet pairs and the branch of mag- 
netic quantum number mp = I atoms. The external field may break a singlet pair into 
two unpaired atoms of mp = 1 under a sufficiently strong magnetic field. A ferromagnetic 
Tonks-Girardeau gas of = 1 atoms appears if the external field exceeds the upper critical 
field Hc2. The singlet pairs and unpaired mp = 1 bosons coexist in the intermediate region 
Hci < H < Hc2. We show that for strong coupling the low energy physics of the gapless 
phase is described by the universality class a two-component Luttinger liquid as long as the 
spin dynamics are frozen out. Moreover, from the TBA equations, we obtain exact results 
for the ground state energy and magnetism for the system with an external magnetic field, 
which provide an exact phase diagram and the universality class of quantum phase transi- 
tions for the integrable spinor F = 1 Bose gas with strongly repulsive and antiferromagnetic 
spin-exchange interactions. 
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II. THE MODEL 



We consider particles confined in ID to a length L with delta-function type density- 
density interactions and spin-spin interactions between two atoms. In first quantized form, 
the Hamiltonian of this model is given by 6| 



^ = E ^ + Et^o + c^F. ■ F,]5(x. - X,) + (1) 

1=1 * i<j 

where Fj are spin-1 operators, cq = {go + 2(72) /3 and C2 = ((72 — fl'o)/3 are the interaction 
parameters which are related to the s-wave scattering lengths ap in the spin-0 and spin- 
2 channels, and gp = Anh^ap/m with m being the mass of the atoms. The term 
accounts for the Zeeman energy which will later be given explicitly. Throughout this paper 
we use the dimensionless units of ^ = 2m = 1 for convenience. We are interested in the 
antiferromagnetic case C2 > and cq = C2 = c where this model is exactly solvable by the BA 
{2^. The model that we examine also has repulsive density-density interactions since the 
interaction parameter cq = C2 > 0. Repulsive interactions result in an effective attraction 
in the two-body scattering matrix for the spin-0 channel and an effective repulsion in the 
scattering matrix for the spin-2 channel. This naturally leads to the formation of singlet 
bound pairs in the spin-0 channel 2^. Due to the existence of the spin exchange interaction 
in the Hamiltonian, the number of particles in a particular spin state (m^ = 1, 0, —1) is no 
longer conserved because spin transmutation is allowed to occur. The scattering between 
two particles of spin irtp = 1 and mp = —1 can produce two particles of spin mp = Q and 
vice- versa. The only conserved quantities are the total particle number N and the total spin 
in the 2;-component 5*^. This model possesses U{1) symmetry for charge conservation and 
SU{2) symmetry which corresponds to spin conservation. For weak interaction and in the 
absence of an external field, the spin dynamics is described by the 0(3) non-linear sigma 
model which can be separated from the BA equations 2l|. The charge sector on the other 
hand is described by collective pairing density fluctuations of free boson fields. 

The BA equations for this Hamiltonian acting on a totally symmetric Bose wavefunction 
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are [20 1 



N M 



exp{ikjL) = -Y\_ ^4:{kj - k) J]^ e^2{kj - A^) 



1=1 a=l 
N M 

l[e2{Aa-ki) = - l[e2{Aa-Ap), (2) 

1=1 13=1 

where j = 1, . . . ,N and a = 1, . . . , M. {kj} is the set of quasimomenta for the particles 
and {Aq,} is the set of spin rapidities that characterize the internal spin degrees of freedom. 
The quantum number M is a conserved quantity satisfying the relation M = N — and 
the function 

en{x) = (3) 
X — ind 

where d = c/4. The energy E = A;| and total momentum p = kj of the system can 
be obtained by solving the coupled BA equations for the sets {kj} and {A^}. 



III. THE TBA EQUATIONS 

In the thermodynamic limit N, L —>■ oo with the ratio N/L finite, the sets of solutions 
{kj} and {Aq,} of the BA equations take certain forms. As mentioned in ref. 20(], the 
kj's and Aq,'s can form complex pairs kj = Xj ± id and Aj = Xj ± id where Xj is real. In 
Fig. 1, we show a schematic configuration of the quasimomenta and spin rapidities for the 
ground state. Notice that each pair of kj's share the same real part as a corresponding pair 
of Aq's. The bound states are associated with a pair of mp = ±1 bosons or two itlf = 
bosons. In the absence of an external field, this bound state is created by the operator 
At = [(aj)2 - 2a[al^]/V3 [1^. In the presence of a sufficiently strong magnetic field, the 



singlet bound state of two rriF = bosons is not energetically favored HI]. In addition to 
that, we also have real kj^s and A-strings of the form A^'^ = A^ + i{n + 1 — 2j)d, j = 1, . . . ,n 
as solutions. The spin-strings characterize the spin wave fiuctuations. In the thermodynamic 



limit, the grand partition function is Z = tr{e~^^^) = e~'^^'^ 23|, l25| where the Gibbs free 
energy G = E + Ez — fiN — TS, chemical potential /i, Zeeman energy Ez = —HS^ and 
entropy S are given in terms of the densities of charge bound states and spin-strings which 
are subject to the BA equations ([2]). 

The equilibrium states are determined by minimizing the Gibbs free energy, which gives 
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rise to a set of coupled nonlinear integral equations - the TBA equations, i.e. 

£^{k) = e -^-H- Tai * In (1 + e-''^^^^^) 
+r[ai - as] * In (1 + e-'^^^^'^) 

oo 

-T J][a„_i + a„+i] * In (1 + e-^-^'^'^) , 



n=l 



62{k) = 2{e - c'2 - /i) + T[ai - as] * In (l + g-^^^"')/^ 

+T[a2 - a4 - ae] * In (l + e'^^^'^)/^) , 
Mk) = nH + r[a„_i + a„+i] * In (l + e-''^'^/^) 

oo 

+T5^T„„*ln(l + e-^"('=)/^). (4) 

m=l 

We present a detailed derivation of the TBA equations in the Appendix. In the above 
equations, the convolution / * g{x) and the functions a„(x) are defined in equations flA12p 



and (lA13p . The function Tmn{x) is also given in the Appendix. The TBA equations are 
expressed in terms of the dressed energies ei{k), e2{k) and 4>n{k) for unpaired states, paired 
states and spin-strings, respectively. The dressed energy ei{k) depends not only on the 
chemical potential /i and the external field H but also on the interactions between unpaired 
bosons and singlet pairs as well as the spin fluctuations characterized by the spin-strings 
(f)n{k)- Physically, the dressed energies measure the energies over the "Fermi surfaces". We 
clearly see that spin fluctuations are coupled to the dressed energy of unpaired mp = 1 
bosons ei{k) through the last term in the first equation in (jl]). The spin fiippings caused 
by thermal fluctuations are described by the last equation where the magnon excitations in 
mp = I bosons are described by an effective ferromagnetic spin-spin interaction. There is 
no such spin fluctuation coupled to the dressed energy of bound pairs due to its spin neutral 
effect. 

In the strong coupling limit, the dressed energies ei{k) and e2{k) marginally depend on 
the pressures of each other, see the dressed energy-dependent terms in the first and second 



21, 



28j. The 



equations. This is similar to the configuration for the attractive Fermi gas 
only difference is that here the unpaired bosons may scatter between themselves whereas in 
the attractive Fermi gas unpaired fermions do not scatter among themselves. The pressure 
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per unit length of the system is derived from the expression p = ~dG/dL as 

rji poo 

271 ' 



oo 



rp /•OO 

+ - / In (1 + e-^^«/^) dfc, (5) 

^ J — oo 

where the first term corresponds to the pressure for unpaired bosons and the second term 
to the pressure for singlet pairs. 



IV. QUANTUM PHASE TRANSITIONS AND MAGNETISM 

Solving the TBA equations (jlj) imposes a formidable challenge due to the involvement of 
infinitely many spin-strings. Here we shall focus on the ground state properties and quantum 
)hase transitions driven by an external magnetic field. Following the method developed in 
26|, |29|], we state two conditions to proceed on, namely we consider: (I) the scenario where 
we are in the ground state with T — > 0, and (II) the strong coupling limit c ^ 1. With these 
two conditions, we can obtain a series expansion in terms of the coupling strength 1/c for 
various thermodynamic quantities, as we shall see later. The strong interaction condition 
should be easily reached because generally the interaction energy is much larger than the 
kinetic energy for a dilute gas in ID with finite interaction strength c. 
When T ^ 0, the TBA equations (jlj) simplify to 



ei{k) = k'^ — 11 — H + a4^* {k) + [a^ — ai] * £2 {k) 
e2{k) = 2(A;2 - c'^ - /i) + [a^ - oi] * £^{k) 

+ [aQ + a^- a2]* e2ik), (6) 

where the dressed energies e~{k) with a = 1, 2 imply that we only consider the domain 
where the function ea{k) < 0. The negative part of the dressed energies ea{k) for k < 
Qa corresponds to occupied states in the dressed energies while the positive part of Sa 
corresponds to unoccupied states. The integration boundaries Qa characterize the "Fermi 
surfaces" at ea{Qa) = 0. There are no A-strings involved in the ground state (all (f)n{k) are not 
occupied), thus the dressed energy equations evolve into two coupled dressed energies. This 
characterizes the scattering among singlet bound pairs and unpaired bosons. They provide 
complete phase diagrams and information of quantum phase transitions with respect to the 
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Zeeman splitting parameter H and the chemical potential /i. The pressure of the system 
can be represented in a neater way if we introduce the following notation when T ^ 0: 

p = / e]^{k)dk / e2{k)dk 

= P1+P2- (7) 

Every thermodynamic quantity with a subscript 1 (or 2) corresponds to unpaired states (or 
paired states). 

When c 1, we can take a Taylor expansion of the functions an- Throughout this paper, 
we only keep track of terms up to order 1/c^. Higher order corrections can be calculated in 
a straightforward manner. In this limit, equation becomes (up to order 1/c^) 

e,{k) ^ k^-^-H-^ + ^-Pl, 
' ^ 2c' 5c'' 

.2(fc)-2(^^-c--^) + |l + ^. (8) 

We then integrate equations ([8]) between the "Fermi points" ±Qi and ±(^2 so that we can 
re- write the equations in terms of pi and p2- This gives 

- 2xp, « \q\ - 2^Q, - 2HQ, + 

-.p. - \qI - ic-Q. - + %2 + ki. (9) 

3 5c' 6c' 

We also make use of the fact that the dressed energies ei{k) and e2{k) vanish at the "Fermi 
points" i.e., £i(±(5i) = and £2(±Q2) = 0, 

Substituting the "Fermi points" into equations ^ and then re-arranging and iterating 
the terms yield 

o 3/2 3/2 / A \2 

2^1 , /^l Pi 4p2 ^ 

Pi ~ -T, \- 



37r 47r \2yUic' 5/iic' 



3/2 / A 
/i/ / pi 4p2 



TT \2yUlC' 5yUlC' 

3/2 3/2 



37r TT \5/i2c' 24/i2c'/ 



4f.r= / + ^ , (11) 



TT \5/i2C' 24/i2C' 



where we denote the effective chemical potentials fii = fi + H and /i2 = /i + c'^ for the 
unpaired and paired bosons. From the relations n = dp/dfi and nm^ = dp/dH and after 
some lengthy iterations, we arrive at the expressions for the chemical potentials of unpaired 
and paired bosons, 



Hi 



2(1 -m')^' 



16m^ 32(1 — 



57 



157 



16 

512fm^)3 



, 4(1 -m") 32m' 



97 



57 



+ 



157 



(12) 



where ■j = c/n. Substituting these two equations back into pi and p2 gives the pressures 

2 



pi ^ -TT n'^(m^)'^ 1 

3 V 7 

P2 ~ — 7r^n^(l - m^)^ I 1 



6m^ ^ 48(1 — m 



57 

(1 - m"^) 48m' 



27 



57 



(13) 



Further, the free energy can be obtained as 
F 



1 2 3/ ZS3 . 4m" 32 1-m" 

-7r^n^(m")^ 1 ^ ^ 

3 ^ ^ V 37 57 

72(1 - m") - Hnm'' + O 



m 



+ 



32m' 
57 



To find the ground state energy, we can use the relation E = F + Hnm^ . There is also 
an alternative way to derive the ground state energy based on the definition E = and 
the distribution of {kj} in quasimomenta space. Indeed, we show that the energy per unit 
length derived from the discrete BA equations for arbitrary magnetization. 



2(32^2 - lOni) 3(32^2 - lOni) 



5c 



+ 



25c2 



1 



+-7r^nl ( 1 + 



2(48ni + 5n2) 3(48ni + Sns 



15c 



+ 



225c2 



(15) 
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coincides with the TBA result E = F + Hnm^ up to the order of 1/c through the relations 
rii = nm^ and n2 = f(l — where rii and n2 are the density of unpaired and paired 



bosons, respectively. However, the dressed energy formalism provides a more elegant way 
to study quantum phase transitions |30 |. 

For strong coupling, a pair of two bosons becomes stable because the binding energy 
~ ""^^ exceed the kinetic energy. Therefore, the ground state in the absence of 

an external field is characterized by an empty "Fermi sea" for unpaired bosons and a fully 
filled "Fermi sea" for bound pairs. From the dressed energy equations ([6]), we find that 
quantum phase transitions driven by an external field can be determined by the energy 
transfer relation H = — fi2 + c^/16, i.e.. 



where we have used the relations /i^ = {E/L + n2eh) for k = 1, 2 to obtain expressions 
for the chemical potentials up to order 1/7^. 

The lower critical field Hd diminishes the gap, thus a phase transition from a singlet 
ground state into a gapless phase, where two dressed energies of the paired and unpaired 
bosons couple to each other, occurs when H > Hd- When the external field exceeds the 
upper critical field Hc2, all singlet bound pairs are broken which leads to a ferromagnetic 
Tonks-Girardeau Bose gas. The lower and upper critical fields are found by letting = 
and = 1 in (|T6l) . with result 



In Fig. 2, we show the magnetization vs external field for different values of the interaction 
strength c. We see clearly that for H < Hd there is no breaking of bound pairs. The 
magnetization gradually increases from zero to n as H gradually approaches Hc2- The 
phase transitions across Hd and Hc2 are of second order. In the vicinities of Hd and Hc2, 




(16) 




(17) 
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the leading order of the respective normahzed magnetizations are given by 

8{H-H^i) / 86 2813 



^ 7r2n2 V 45O72 

(H^2 -H) f 72 2536 \ 

ml ^ I - ^ , , 1 + ^ . 18) 

2 27r%2 \^ ^ 57 2572 ) ^ ' 

which show a hnear dependence on the external field near the critical points. For an external 
field Hci < H < Hc2, the singlet paired state and unpaired state coexist. They form a two- 
component Luttinger liquid in this gapless phase. 

In Fig. 3, we show the ground-state phase diagram in the n — H plane. As n — > 0, both 
critical fields approach the same value = eh/2. The solid (dashed) lines correspond to 
the two critical fields for the case c = 20 (c = 40). The ferromagnetic phase of all atoms in 
state \F = l,mp = 1) appears above the critical field Hc2, the singlet phase of singlet pairs 
appears below the critical field Hd and the mixed phase of atoms in state \F = l,mF = 1) 
and singlet pairs appears between the two critical fields. 

V. THE SPIN AND CHARGE VELOCITIES 



In ID systems, spin-charge separation is the hallmark of many-body physics 3l|]. The 



collective charge excitations are described by sound modes with a linear dispersion. The spin 



excitations are gapped with a dispersion e^{p) = A/A2~+^^2p2 where is the excitation gap 
and v,y is the spin velocity in spin branch u. This leads to the phenomenon of spin-charge 
separation. A method has been proposed to probe this phenomenon experimentally in a ID 



system of interacting electrons at low energies [32 1. 

To calculate the charge velocity, we need to find the energy of the lowest excited state that 
does not involve breaking any pairs. In the absolute ground state where H = 0, the system 
is only made up of fully paired states below the "Fermi level" and the total momentum of 
the system is zero. This is achieved when there is no magnetic field present. To excite the 
system, we allow the pair with the largest momentum to leave the "Fermi sea" and let the 
excited state have a total momentum of p, i.e., Y2j = P- We then calculate its total energy 
E = k'j. The difference between the excitation energy and the ground state energy is 
equal to the charge velocity times p. The energy difference in the thermodynamic limit is 
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thus calculated to be 

E{p) -Eo = nn,p (^1 + ^^+0 (^i^ . (19) 

Therefore in terms of the total particle number n = N/L and interaction strength c, the 
charge velocity is 

An alternative way to calculate the charge velocity for the singlet ground state is through 
the relation 




(21) 

y 74 y uij' J 

Both methods yield the same result. 

The spin velocity on the other hand is calculated by considering the lowest excited state 
where one pair is broken into two unpaired states. Both unpaired states will occupy opposite 
ends of the momentum distribution so that the excitation energy is minimized. The total 
momentum of the excited state can be parameterized by p in the same manner as before. 
We can equate the energy difference between the excited state and the fully paired ground 
state to the energy dispersion e(p). In the thermodynamic limit, 



(22) 



From the original relation e(p) = a/A^ + vl^P^^ we obtain the relation e(p) = A + in the 
limit A ^ 1 where the gap is very large. Comparing both expressions for the dispersion 
energy, we can easily verify that A = c^/8 and 

-2 ^ 64^2 \ 



5c 



- ) ■ (23) 



Hence the spin velocity is 



The spin velocity is divergent due to a large energy gap as c — oo. This demonstrates 
that there is spin-charge separation over the singlet ground state. We also note that this 
phenomenon depends on the state of the system within an external field. Essler et al. 21 1 
showed that spin and charge velocities are equal in the weak coupling limit when there is 
no external field involved. 
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In the gapless phase when Hd < H < Hc2, the ground state (T — > 0) of this system is 



conformally invariant 



33 



3J] . The excitations close to the "Fermi surfaces" in unpaired and 
pair branches have hnear dispersions. The finite-size corrections to the ground state energy 
are given by 

Eo = Le^ - ^{vi'^ + vi'^) (25) 

where the central charge c = 1 for this system, Eq is the ground state energy for the 
finite system and is the ground state energy density for the infinite system. The charge 
velocities for unpaired and paired bosons are given explicitly by the expressions 

(1) ^ / , 2(32^2 -lOm) , 3(32n2 - 10ni)2 

(2) A , 2(48ni + 5n2) , 3(48ni + 5n2)2\ 

In this phase, spin fluctuations are frozen out and thus the charge density fluctuations 
dominate the ground state and is effectively described by the universality class of a two 
component Luttinger liquid. 



VI. CONCLUSION 

In conclusion we derived the TBA equations for a system of ID spin-1 bosons with 
repulsive density-density and antiferromagnetic spin exchange interactions and solved the 
TBA equations for the zero temperature case in the strong coupling limit. We obtained 
the ground state energy, chemical potentials, critical fields and magnetization in terms of 
interaction strength and the external magnetic field. We also presented an exact phase 
diagram of strongly interacting spin-1 bosons which facilitates experimental analysis of phase 
segments. For the weak coupling limit, the collective excitations in the charge sector is 
described by a Tomonaga-Luttinger liquid, whereas the spin dynamics is described by the 



0(3) non-linear sigma model 2l|. However, for the strong coupling limit, spin fluctuations 
can be suppressed by a strong external magnetic field. The density fluctuations thus evolve 
into a two-component Luttinger liquid. At zero temperature, the model exhibits three 
quantum phases: singlet pairs of two bosons for external field H < Hd] a fully-polarized 
Tonks-Girardeau gas phase of mp = 1 bosons for H > Hc2] and a mixed phase of singlet 
pairs and unpaired mp = 1 atoms for an intermediate field Hd < H < Hc2- The phase 
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transitions in the vicinities of Hd and Hc2 are of second order with a hnear-field- dependent 
magnetization. Our results provide a new aspect of this model, namely spin liquid v.s. 
Luttinger liquid behavior. 



This work has been supported by the Australian Research Council. We thank Profs J.-P. 
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APPENDIX A: DERIVATION OF THE TBA EQUATIONS 



Here we derive the TBA equations in detail following the steps for the attractive spin-1/2 



fermion model in Chapter 13 of Takahashi's book 



25| . Substituting all possible solutions 



for {kj} and {Aq,} back into the BA equations ([2]) gives 

i^j 1=1 

00 M„ 



n=l a=l 
00 Mn 



(Al) 



n=l a=l 



exp 



'i2\jL) = JJe_i(Aj— A;i)e5(Aj - A;, 



1=1 
N2 



]^e_2(Aj - Az)e4(Aj - Xi)eQ{Xj - A; 



(A2) 



where Emn{x) is defined as 



J]e(„_i)(A^-A;Oe(„+i)(A^-A;z) 
1=1 

00 Mm 

m=l (3=1 



e|m,-n|(a;)e^^_„l^2(^) • • ■ 
el(x)el{x) . . . 



'x)e2n{x), 



for n = m. 



(A3) 



(A4) 
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Following the technique pioneered by Yang and Yang 23|] for spinless bosons, the loga- 
rithm of equations (Al) to (A3) then gives 



1=1 

00 



1=1 ^ 

00 Mn 



kj- \i\ ^ ( kj - Xi 



5c' 



n=l a=l '- 



2XjL = 277 Jj + 



1=1 



{n - l)c' 
Xj - ki 



n + l)c' 
Xj - h 



5c' 



N2 



A, — A; 



2c' 



A, - A^ 
4c' 



6c' 



E 



A!!-fc 



[n - l)c' 

00 A/i 



27rL^^ + ^5^e„„ 

m=l /3=1 



A^-fc/ 
(ri + l)c' 



\ n Am 



where ^^fx) = 2 tan ^ix) and 



+ 2^ 



+ 



©mn (•^) ^ 



(A5) 



(A6) 



(A7) 



(A8) 



|m— n| j \ |m— n|+2 

2^^ (f ) + 2^^ (f ) + . . . 
,+2^fe)+^(l^)' forn = m. 

Writing the occupied distribution functions of the unpaired fc's, the paired fc's and the A- 
strings as Pi{k), P2{k) and cr„(A;) and their corresponding unoccupied distribution functions 
as Pi{k), P2{k) and we take the thermodynamic hmit of the above equations to obtain 

the integral equations 

1 

2^ 

+ [05 - ai] * p2{k) 



pi{k)+p\{k) = — + ai*piik) 



^[a„_i + a„+i] * anik), 



(A9) 



n=l 
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P2{k)+p^{k) = - + [og - ai] *pi(A;) 



TT 



+ [a6 + a4 - 02] * P2(A:), 



(AlO) 



an{k) + (7„(A;) = [a„_i + an+i] * pi(A;) 

00 



m=l 



where 



f{x — x')g{x')dx' , 



n\c 



TT (nc')^ + X 



2 ' 



(All) 

(A12) 
(A13) 



and 



-^mn('^) 



a|„_„l(x) + 2aim-n\+2ix) + 

-\-2am+n~2{x) + am+n(2^)) 

202(2;) + 2a4(x) + . . . 
^ +2a2n-2(a;) + a2n{x), 
The distribution functions are related to the particle numbers via the relations 



for n ^ m, 



for n = m. 



(A14) 



712 = 



L 
N2 



m 



L 

00 

= 1: 



CX) 

n / 

n=l 



Pi{k)dk, 
P2{k)dk, 

an{k)dk. 



(A15) 
(A16) 

(A17) 



The total number of microstates in an interval dk is 

{L{p,{k) + p^l{k))dk)\ 



dW 



{Lpi{k)dk)\{Lp'l{k)dk)\ 
{L{p2{k) + pm)dk)\ 
{Lp2{k)dk)\{Lp^{k)dk)\ 

{L{an{k) + a^{k))dk)\ 



X 



n 



\ {Lan{k)dk)\{La^{k)dk)\' 



(A18) 
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Through the use of Stirhng's approximation, the entropy is written as 

[(pi + p\) ln(pi + pt) 

) 

-pi Inpi - p^ Inp^] dk 



S 
L 



/oo 
[(P2 + p^) ln(p2 + p^; 
-oo 

-P2 In p2 - p2 In P2] (ifc 

n=l 

-a„ In an - cr'^ In c?A;. (A19) 
The Gibbs free energy per unit length is 

n = ^-pn-^ + ^ (A20) 
where /x is the chemical potential, the energy per unit length is 

771 noo POO 

= epi{k)dk+ {2k^ -2c'^)p2{k)dk (A21) 

^ J —00 J —00 

and the Zeeman energy per unit length is given in terms of the external magnetic field H 

rp poo °° POO 

-y = -H j pi{k)dk + H^nj an{k)dk. (A22) 
i-00 J-00 



Minimizing the Gibbs free energy and going through a similar procedure as shown in 25| . 
we arrive finally at the TBA equations (jl]). 
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FIG. 1: (color online) Schematic configuration of the quasimomenta k and spin rapidities A in 
the complex plane for the ground state in the strong coupling limit with = 17 and M = 12. 
Dashed boundaries indicate pair formations with the same real parts. Each pair in fc-space has 
a corresponding pair in A-space. As mentioned in the text, the individual particle numbers for 
mp = are not conserved and thus can fluctuate. 
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FIG. 2: Normalized magnetization versus re-scaled magnetic field H/eb for different values of 
interaction strength c. 




FIG. 3: Phase diagram in plane (n, H) for different values of interaction strength c. 



21 



